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ABSTRACT 

We present hydrodynamical iV-body simulations of clusters of galaxies with feedback taken 
from semi-analytic models of galaxy formation. The advantage of this technique is that the source 
of feedback in our simulations is a population of galaxies that closely resembles that found in the 
real universe. We demonstrate that, to achieve the high entropy levels found in clusters, active 
galactic nuclei must inject a large fraction of their energy into the intergalactic/intracluster media 
th roughout the growt h period of the central black hole. These simulations reinforce the argument 



of Bower et al. (2008), who arrived at the same conclusion on the basis of purely semi-analytic 



reasoning. 

Subject headings: hydrodynamics - methods: N-body simulations 
ies:) cooling flows - X-rays: galaxies: clusters. 



galaxies: clusters: general - (galax- 



1. INTRODUCTION 

X-ray observations of clusters of galaxies al- 
low us to probe the physical properties of the 
hot, diffuse intracluster medium (ICM). In the 
simplest self-similar collapse scenario, where the 
ICM is heated solely by gravitational processes 
(such as adiabatic compression and shocks in- 
duced by supersonic accretion), we expect the 
X-ray luminosity Lx of clusters to scale with 
gas tempe rature T as L x °c T 2 (for T > 2 
keV); see iKaiserl (Il986l) . However, the ob- 
served luminosity-temperature relation is much 
steeper than predicted, with Lx oc T 2 5 ~ 3 at 
T > 2 keV (lMarkevitcb|[l99l lArnaud fc Evrard 
19991 : IWu et all Il999l : lEttori et all l2002h . be- 
comi ng steeper still at group scales T < 2 



due to an exc e ss of e ntropy [j] in cluster cores 



( Ponman etal] Il999t iLlovd-Davies et al.l [2000; 



Finoguenov et al.l 120021 ) . The source of this ex- 
cess entropy is likely to be a combination of non- 
gravitational cooling and heating processes (see 



Voit 2005 for a review). 

Radiative cooling leads to an increase in the 
entropy of the ICM since cold, low-entropy gas 
condenses to form stars, reducing the gas density 
and causing higher-entropy gas to flow inwards to 
compensate for the loss of pressure support. This 
leads to a decrease in X-ray luminosity, breaking 
the self-similarity of the cluster scal ing relations 



comi ng steeper still at group scales 1 >, z 
keV dHelsdon fe Ponmanl l200d IXue fc Wul l200d: 
lOsmond fc Ponmanl 12004 ) . This deficit in X-ray 
luminosity, particularly in low-mass systems, is 



Pearce et al. 


20001 iMuanwone et al. 2001, l2002t 


Dave et al. 


2002t IVoit et all 20021 Wu & Xue 



2002). However, hydrodynamical TV-body sim- 



ulations that only include radiative cooling gener- 



1 We define entropy as S = kT /nZ , where k is Boltzmann's 
constant, n c is the electron number density and 7 = 5/3 is 
the ratio of specific heats for a monoatomic ideal gas. 
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ically predict an exc essive amount of cold gas that 
can form stars (e . g. ISuginohara fe Ostrikerlll998t 
iLewis et all l2000t lYoshida et al.l l2002h : typically 
the fraction, /*, of baryons in stars lies in the range 
0.3 — 0.5 in simulated clusters, in clear conflict 



with t he observed val u e of f* < 0.1 ( Balogh et al 



20011 : iLin et alJ 120031 : Balogh et al 



20081 ). This 



highlights the need for extra heating (feedback) 
from astrophysical sources to regulate the cooling 
flow and quench star formation. 

Non-gravitational heating, occurring before or 
during gravitational collapse, raises the entropy of 
the ICM, preventing gas from reaching high den- 
sities in central cluster regions and reducing its X- 
ray emissivity. This effect will be greater in lower- 
mass systems, leading to a steepening of the Lx-T 
relation as desired. The most obvious sources of 
non-gravitational heating are Type II supernovae 
(SNe) and active galactic nuclei (AGNs). In the 
simplest heating models, the energy released by 
these phenomena is assumed to be injected impul- 
sively into the gas at high rcdshift. Although not 
well motivated physically, these so-called preheat- 
ing models have proved capable of reproducing 
the observed slope and normalisation of cluster 
scaling relations when inco r porated in simulations 
(e.g. iBialek et al.l 1200 it 



Brighenti 



Mathewsl 

20011: iMuanwong et al l 120021: IBorgani et all 1200? 
Tornatore et al.ll2003t IBorgani et alJl2005h . How- 



ever, the preheating scenario suffers from sev- 
eral problems. For example, the predicted scatter 
about the mean Lx-T relation is much smaller 
than observed and large isentropic cores are gen- 
erated in low-mass systems that are n ot seen (e.g. 
Ponman et aDl2003tIp"ratt et alJl2006h . 



Recently, attention has shifted to more realistic 
models which attempt to couple cooling, star for- 
mation and black hole growth with feedback from 
SNe and AGNs. There has been considerable ef- 
fort to include such processes in hydrodynamical 
simulations of galaxy groups and clusters. How- 
ever, an explicit treatment is unfeasible since star 
formation, black hole growth and associated feed- 
back all occur on scales much smaller than can 
be resolved with present computational resources. 
The only option is to develop phenomenological 
prescriptions and assess their validity by compar- 
ing the properties of simulated clusters with as 
much observational data as possible. 

Various theoretical models of stellar feedback 



have been proposed and can be grouped in two 
basic categories: thermal, where the available 
supernova energy is used to raise the tempera- 
ture of neighbouring gas particles, and kinetic, 
where neighbouring particles are given a veloc- 
ity 'kick'. Cosmological simulations with radia- 
tive cooling, star formation and supernova feed- 
back yield ICM profiles and scaling relations that 



both thermal (e.g. Valdarninil 2003; Kav 2004: 


Kav et al. 


2004. 2007; Nagai et al. 


2007) and ki- 


netic (e.g. 


Borgani et al.l 


2004, 2005 


; Romeo et al. 



2006) models. In particular, the mean Lx-T rela- 



tion is generally well reproduced on cluster scales 
T > 2 keV. However, X-ray luminosities are found 
to be substantially larger than the observed values 
on group scales. Other problems are that observed 
baryon fractio ns are typically smaller than in sim- 



ulatio ns (e.g. iKravtsov et alJ 120051 : lEttori et al 



2006), and star formation is still too efficient, in- 
dicating that additional feedback mechanisms are 
required to offset radiative cooling. 

The favoured candidate is the gravitational 
energy liberated by the accretion of gas onto 
central supermassive black holes within galax- 
ies. This can be extremely large, of the order 
1Q 62 (Mbh/1Q 9 M q ) erg, where M BH is the black 
hole mass. Analytical calculations suggest that 
the Lx-T relation of groups can be accounted for 
if ~ 1% of the energy released b y AGNs is cou- 



pled to the surrounding ICM (e.g. ICavaliere et al 



2002). The precise details of how this energy is 
transferred to the ICM are not well understood 
at present, but it appears there are two major 
channels via whic h black holes interact with their 
surroundings (see lMcNamara fe Nulsenll2007l for a 
review). 

At high redshift, mergers of gas-rich galax- 
ies occur frequently and are expected to fun- 
nel copious amounts of cold gas towards galac- 
tic centres, leading to high black hole accretion 
rates and radiating enough energy to support 
the luminosities of powerful quasars. Quasar- 
induced outflows have been obser vationally con 



firmed in a number of cases 



e^ 



Chart as et al 



2003t ICrenshaw et~aTl l2003t fPounds et alJ 120031 ). 



and demonstrated in high-resolution simulations 
of galaxy mergers. In such simulations it is usually 
assumed that a small fraction of the bolometric lu- 
minosity can couple thermally to the surrounding 
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gas. This approach has been shown to success- 
fully repr oduce the observed black hole-bulge mass 



relatio n (|Di Matteo et al.l 120051 : iRobcrtson et al 



2006a) and explain the red colours (ISpringel et al 



iRober 
s (Spri 



2005a[ ), X-ray haloes (ICox et all 120061) an d fun 



damental plane relation (IRobertson et al.l 120 06b) 



of ma ssive elliptical galaxies. iDi Matteo et al 
(|2008l ) recently extended this work to fully cos- 
mological simulations, obtaining similarly encour- 
aging results. 

There is another mode of AGN feedback which 
is not related to quasar activity or triggered by 
galaxy mergers. Evidence for this can be seen, 
for example, in brightest cluster galaxies (BCGs) 
which contain X-ray cavities filled with relativis- 



tic plasma (e.g. Blanton et al. 200ll Birzan et alj 



2004t iMcNamara et al.ll2005t iFabian et al.ll2006f h 

These radio-loud X-ray depressions, referred to 
as bubbles, are thought to be inflated by rela- 
tivistic jets launched from the central black hole. 
Bubbles may rise buoyantly, removing some of 
the central cool gas and allowing it to mix with 
hotter gas in the outer regions of groups and 
clusters. Together with the accompanying me- 
chanical heating, this can constitute an efficient 
mechanism for suppressing cooling flows. Sim- 
ulations of idealised clusters, performed with 
hydrodynamical mesh cod es, suggest that this 
is indeed the case (e.g. IChurazov et al.l 12001 



Quilis et al. 2001; Ruszkowski & Begclman 2002; 
Bruggen et alj2002tlBruggenl2003l ; lDalla Vecchia et al 



20041) . 

The first attempt to implement a self-consistent 
model of black hole growth and heating by AGN- 
driven b ubbles in cosmolog i cal si mulations was 
made by ISiiacki fc Springell (|2006l ). They found 
that bubble injection can substantially affect the 
properties of the ICM, especially in massive, re- 
laxed systems at late times. In particular, bub- 
bles are able to efficiently heat central cluster gas, 
reducing the amount of cold baryons and star for- 
mation in the central cD galaxy. Furthermore, the 
gas density is reduced and the temperature is in- 
creased out to radii ~ SOOh^ 1 kpc, leading to a 
decline in X-ray luminosity and establishing a flat 
entropy profile in central regions. Such trends are 
precisely what is required to reconcile simulations 
of galaxy clusters with observations. 



corporating both 'quasar' and 'radio' modes, as 
well as star formation and feedback from SNe. Mo- 
tivated by the properties of X-ray b inaries (e.g. 
Fender et all 1 19991: iGallo et alj 120031) . the tran- 



sition between the two states is assumed to be 
governed by an accretion rate threshold Xradio = 
Mna/MEdd, where MEdd is the Eddington ac- 
cretion rate. Quasar-like feedback occurs for ac- 
cretion rates greater than Xradio, while mechan- 
ical bubble heating takes place otherwise. This 
model has recently been used by IPuchwein et al" 



Following this work, ISijacki et alJ (|2007l ) have 
formulated a unified model of AGN feedback, in- 



(|2008l ) to investigate the Lx-T relation and gas 
fractions of galaxy groups and clusters in cos- 
mological simulations. They demonstrated that 
AGN feedback reduces the X-ray luminosities of 
groups and poor clusters more than rich clusters, 
resulting in a steepening of the Lx-T relation on 
group scales. In fact, the X-ray properties of their 
simulated objects are in excellent agreement with 
observational data on all mass scales. However, 
since their sample size is quite small (21 objects), 
it is unclear whether their model can generate a 
realistic population of cool core (CC) and non- 
cool core (NCC) systems and thus explain the 
observed scatter about the mean Lx-T relation. 
The gas fraction in groups and poor clusters was 
shown to decrease significantly with the inclusion 
of AGN heating (even though fewer baryons were 
converted into stars), because gas is driven from 
central regions to cluster outskirts. In more mas- 
sive systems, the main effect of AGN feedback is 
to lower the central gas density and substantially 
reduce the number of stars formed. However, even 
with stellar and AGN feedback, the stellar fraction 
within the virial radii of their simulated objects 
appears larger than observations suggest, imply- 
ing that the cooling flow problem has not been 
fully resolved. 

In this paper we pursue a different, but com- 
plementary, approach to the theoretical study of 
groups and clusters of galaxies. Instead of under- 
taking fully self-consistent hydrodynamical simu- 
lations, we investigate what current semi-analytic 
models (SAMs) of galaxy formation predict for the 
thermodynamical properties of the ICM. Our goal 
is to extend the predictive power of these models, 
thus providing additional constraints on, and in- 
sights into, the physics of galaxy formation. We 
hope that our work will be useful for guiding the 
development of future SAMs that can simultane- 
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ously account for the properties of both the galaxy 
distribution and the ICM. 

The basis of a SAM is a set of dark matter 
halo merger trees, usually obtained from a high- 
resolution TV-body simulation. The behaviour of 
baryonic matter within these dark haloes is mod- 
elled using analytical 'recipes' to capture the es- 
sential physical processes involved in the forma- 
tion and evolution of galaxies: gas cooling, star 
formation, black hole growth, feedback, galaxy 
dynamics, galaxy mergers, etc. These recipes typ- 
ically contain a number of adjustable parame- 
ters, which are tuned to attain the best possi- 
ble match to selected observational properties of 
the galaxy distribution. This semi-analytic ap- 
proach has proved largely successful, reproducing 
many key properties of real galaxies such as lumi- 
nosities, colours, star formation rates, the Tully- 

Fisher relation and the black hole-bulge mass 

I 1 1 — ' — I I — 1 

relation (e.g. Kauffmann et alj 119 99; Cole et al 

20001: iKauffmann fc Haehneltl l2000l; iBenson et al 



2003; Dc Lucia et al. 2004; Granato et al. 2004; 



Bower et al.ll2006l:ICroton et al.ll200dlCattaneo et alJ 
2006[lMenci et al.l2006tlMonaco et al.l2007tlFont et alJ 

2008). 



The essence of our method is to generate a 
semi-analytic galaxy catalogue for a cosmologi- 
cal iV-body simulation, compute the energy re- 
leased by SNe and the AGN in each model galaxy, 
and inject this energy at the appropriate position 
and time into the baryonic component of a non- 
radiative hydrodynamical simulation that has the 
same distribution of dissipationless dark matter. 
In this way we can track the effect of energy feed- 
back from semi-analytic galaxies on the intraclus- 
ter gas. A similar technique has already bee n used 
to study the metal enr ichment of the ICM (jCora 
200dlCoraet~aTll2008h . 

There are several benefits of our hybrid ap- 
proach. Firstly, feedback is guaranteed to orig- 
inate from a galaxy population whose observa- 
tional properties agree well with those of real 
galaxies. This is generally not the case in self- 
consistent hydrodynamical simulations. Secondly, 
we only need a single dark matter simulation to 
construct a semi-analytic galaxy catalogue. Since 
A^-body simulations are much less computation- 
ally demanding than hydrodynamical simulations, 
we can, in principle, attain significantly greater 
resolution for the collisionless component. With 



a high-resolution dark matter simulation we can 
construct a comprehensive galaxy catalogue and 
thus a detailed model for feedback from galax- 
ies. Thirdly, the energy transferred to the ICM 
by SNe and AGNs can be calculated directly from 
the semi-analytic galaxy catalogue, before cou- 
pling the SAM to a hydrodynamical simulation. 
By avoiding the need to include an explicit sub- 
grid model for star formation, black hole growth 
and associated feedback processes, the hydrody- 
namical simulations we use to track the injection 
of energy into the ICM require considerably less 
computational effort. Finally, as we shall see, a 
lower resolution can be used for the gas than the 
dark matter in our hydrodynamical simulations. 
As a result, our technique is readily applicable to 
large cosmological volumes, allowing massive clus- 
ters to be simulated. 

The layout of this paper is as follows. In Sec- 
tion [2] we discuss our numerical method in detail, 
describing the simulations used, the relevant com- 
ponents of the SAM and our feedback implemen- 
tation. We use our hybrid approach to investigate 
the bulk properties of the ICM in Section [3l com- 
paring our results with a selection of observational 
data. A resolution test is presented in Section [|] 
to demonstrate the robustness of our results. We 
summarise our results and conclude in Section [5] 

2. THE NUMERICAL MODEL 

Our hybrid technique for studying the effect 
of galaxy feedback on intracluster gas consists of 
three distinct components: an underlying dark 
matter simulation, a semi-analytic galaxy cata- 
logue built on the halo merger trees of this simu- 
lation, and a hydrodynamical simulation to track 
the energy injection from model galaxies. In this 
section we provide a detailed description of each 
part of the modelling process. 

2.1. Dark matter simulations 

The cosmological model we adopt is a spatially- 
flat ACDM model with cosmological parameters 
fi m ,o = 0.25, n Afi = 0.75, h = 0.73, n s = 1 
and as. = 0.9. Here fl m ,o and ^a,o ar e the to- 
tal matter and dark energy density parameters, h 
is the Hubble parameter Hq in units of 100 km 
s _1 Mpc -1 , n s is the spectral index of primor- 
dial density perturbations, and <7s,o is the rms 
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linear density fluctuation within a sphere of ra- 
dius 8/i _1 Mpc. The subscript signifies the 
value of a quantity at the present day. While 
there is some tension between our chosen param- 
eter values (particularly n s and as,o) and those 
derived from the five-year Wilkins on Microwave 



Anisotropy Probe (WMAP) data (jPunklev et al 



2009). we have deliberately used the same cosmo- 
logical parameters as the Millennium simulation 
(jSpringel et al. 2005bh since we eventually hope 
to apply our technique to the full Millennium 
volume, using the publicly-avai lable Millen nium 
semi- analytic galaxy catalogues ( Lemsonll2006l ) as 
input. 

Initial conditions were created at a redshift z\ — 
127 by displacing particles from a glass-like distri- 
bution, so as to form a random realisation of a 
density field with a ACDM l inear power spectrum 
obtai ned from CMBFAST (jSeliak fc Zaldarriaga 
1996). The amplitudes and phases of the initial 
density perturbations were chosen to be the same 
as those of the Millennium simulation initial condi- 
tions. We generated initial conditions for two cu- 
bic simulation volumes: one with a comoving side 
length L = 62.5h~ 1 Mpc and A DM = 270 3 dark 
matter particles, and the other with L = 125/i _1 
Mpc and -/Vdm = 540 3 . The mass of a dark matter 
particle is then tu-dm — 8.61 x 10 8 h~ 1 M & , as in 
the Millennium simulation. The smaller of these 
simulations is for the purposes of initial model 
discrimination, while the larger simulation is to 
demonstrate our best model in a more cosmologi- 
cally interesting volume. 

The massi vely paral l el tre e iV-body/SPH code 
GADGET-2 (|Springell l2005h was then used to 
evolve the initial conditions to z = 0, with full par- 
ticle data stored at the 64 output redshifts of the 
Millennium simulation: z@3 — 127, zq2 = 80, z$i = 
50, z 6a = 30 and log 10 (l + z n ) = n(n + 35)/4200, 
< n < 60. The Plummer- equivalent gravita- 
tional softening length was fixed at e = 40/i -1 kpc 
in comoving coordinates until z = 3, then fixed in 
physical coordinates thereafter. Note that gravi- 
tational forces were softened on a fixed comoving 
scale of e — 5/i _1 kpc in the Millennium simula- 
tion, a factor of two smaller than our z ~ soft- 
ening length. We chose a larger softening scale for 
two reasons. Firstly, by comparing the number 
of dark matter substructures (see below) formed 
in simulations with e = hhr 1 kpc and e = 10/i -1 



kpc, we found that the more aggressive softening 
scheme yields fewer low-mass objects because of 
two-body heating effects. This means we can con- 
struct a more detailed semi-analytic galaxy cat- 
alogue by setting e = 10ft- -1 kpc. Secondly, we 
will be incorporating gas particles into our dark 
matter simulations and it has been shown that 
an optimal choice for the softening length in hy- 
drodynamical simulations is approximately 4% of 
the mean inter-particle spa cing, corresponding to 
e w lOfe " 1 kpc in our case (jThomas fe Couchman 



199a Borgani et al.ll2006h . 



2.1.1. Dark matter haloes, substructure and 
merger tree construction 

Dark matter haloes are identified as virialised 
particle groups within the simulations using the 
friends-of-friends (FOF) algorithm. We adopt a 
standard FOF linkin g length of 20% o f the mean 
particle separation ( Davis et al. 1985f) and only 
save groups that contain at least 20 particles, so 
that the minimum halo mass is 1.7 x \Q w h~ l MQ. 
FOF group catalogues are produced on the fly and 
in parallel by the simulation code. An improved 



versio n of the SUBFIND algorithm ( Springel et al 



20011 ) is then applied in post-processing to these 
group catalogues to find gravitationally-bound 
dark matter substructures orbiting within the 
FOF haloes. 

To compute a virial mass estimate for each 
FOF halo, a sphere is grown about the min- 
imum of the gravitational potential within the 
group until the mean overdensity enclosed reaches 
Ap cr (z), where A is the desired density contrast, 
Pcr(z) = 3H 2 E(z) 2 /8nG is the critical density and 
E(z) 2 = f7 m ,o(l + z) 3 + ^A,o in a ACDM cosmolog- 
ical model. The mass enclosed within the sphere, 
Ma, is then related to the sphere radius, ta, and 
the circular velocity, «a, at this radius by 



Air 

M A = yr A Ap cr (z) 



47rG 3 Ap cr (z) 



1/2 



V A . 



(1) 

Following Croton et al. ( 20061 ) , a density contrast 
A = 200 is used to define the virial mass, radius 
and velocity of a halo throughout this paper. 

Once all haloes and subhaloes have been iden- 
tified for each simulation output, we construct 
merger trees that describe how haloes grow as 
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the universe evolves. This is done by exploiting 
the fact that each halo will have a unique descen- 
dant in a hierarchical s c enario of structure forma- 
tion; see lSpringel et al.l (|2005bl ) for further details. 



Each individual merger tree then contains the full 
formation history of a given halo at z = 0. 

2.2. The semi-analytic model of galaxy 
formation 

Dark matter halo merger trees form the back- 
bone of modern hierarchical models of galaxy for- 
mation. We have generated galaxy catalogues 
for our two dark matter simulations by apply- 
ing a SAM to to the merger trees. The SAM 
we use is the highly s uccessful Munich L-Galaxie s 
model described by iDe Lucia fc Blaizot ( 2007 ). 
We adopt the same set o f model parameters as 
De Lucia fc Blaizotl (|2007h since these parameters 
were used to produce the publicly-available Millen- 
nium galaxy catalogue. As in the full Millennium 
catalogue, galaxy and host halo properties are 
stored at the same 64 redshift values as the sim- 
ulation outputs. We now review the components 
of this model that are relevant for our work; for a 
ful l description o f L-Ga laxies we refer the reader 
to Croton et al. (|2006l) and IDe Lucia &: Blaizotl 



J2007I) 



2.2.1. Star formation and supernova feedback 

The infall of gas into the potential well of a dark 
halo causes the gas to be shock heated to the halo 
virial temperature. At late times and in massive 
systems, this gas is added to a quasi-static hot 
atmosphere that extends roughly to the virial ra- 
dius of the dark halo. Gas from the inner regions 
of this atmosphere can then cool and accrete onto 
a central cold gas disk. At early times and in lower 
mass systems, the post-shock cooling time is suffi- 
ciently short that a quasi-static halo cannot form. 
Instead the shocked gas rapidly cools and settles 
onto the cold disk. 

Once gas in the cold disk exceeds a critical 
surface density, it can collapse and form stars 
(|KennicuttJll989l ). Massive stars rapidly complete 
their life cycle and explode as Type II SNe, in- 
jecting gas, metals and energy into the surround- 
ing medium, reheating cold disk gas and possibly 
ejecting gas from the quasi-static hot halo. For 
a given mass Am* of stars formed over some fi- 



nite time interval, the amount of energy released 
by Type II SNe that is available for heating gas is 
approximated by 



A-Esn = -ehaioAm*w| N , 



(2) 



where «sn = 630 km s^ 1 , based on a standard 
stellar initial mass function (IMF), and ehaio is 
an efficiency parameter. The published value fo r 
this parameter is ehaio = 0.35 (|Croton et al.ll2006h . 
We note that the mass actually locked-up in stars 
is given by AM* = (1 — R)Am*, the rest is as- 
sumed to be instantaneously returned to the cold 
disk. Here R is the recycle fraction, which is as- 
signed the value R = 0.43 in ac cordance with 
the Chabrier IMF dchabrierl l2003h employed by 
De Lucia fc Blaizotl (|2007h . 

The mass of cold gas reheated by SNe is mod- 
elled as 



Am re heated = EdiskAm,, (3) 

where edisk is a parameter that is set to enisk = 
3.5, m otivated by the observational work of lMartin 
(1999). ^f the reheated gas were added to the hot 
halo without changing its specific energy, its total 
thermal energy would increase by 



Ai?hot — 2^ TOrchcatcdU vir- 



(4) 



If there is any excess energy after reheating, i.e. 
AEsn > AShot, then it is assumed that a mass 



A?71 c j ec t e d 



fhalo 



VSN 
V v ir 



£disk 



Am* 



(5) 



of hot gas is ejected from the halo into an external 
'reservoir'. When Ai?sN < Ai^hot, there is insuf- 
ficient energy to eject any gas out of the halo and 
Am c j cct ed is set to zero. Note that the reheated 
mass is not reduced in this case. It follows from 
equation ((5|) that no hot gas can be expelled if 
v V ir > (ehaio/edisk) 1/2 «SN ~ 200 km s~\ whereas 
the entire hot halo can be ejected for small v v i T . 

2.2.2. Black hole growth and cooling flow sup- 
pression 

The growth of supermassive black holes in L- 
Galaxies is driven primarily by quasar mode ac- 
cretion, as a result of galaxy mergers. Black holes 
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can grow either by merging with each other or by 
the accretion of cold disk gas. The coalescence of 
black holes is modelled simply by taking the sum 
of the progenitor black hole masses. The gas mass 
AMbh,Q accreted during a merger of galaxies with 
respective masses M sat and M ccntra i is 



AM 



BH.Q 



/BH (M sat /M CC ntral)M co ld 

1 + (280 km s^V^vir) 2 ' 



(6) 



where M co id is the total cold gas mass present and 
the constant /bh = 0.03 is chosen to reproduce 
the observed local black hole-bulge mass relation. 
Black hole accretion occurs during both minor 
(M sat < M contra i) and major (M sat « M contra i) 
mergers, although the efficiency in the former case 
is reduced by the M sat /M contra i term. 

Once a static hot halo has formed around the 
host galaxy of a black hole, it is assumed that there 
is also continual and quiescent accretion onto the 
central black hole directly from the hot phase. The 
growth rate Mbh.r of the black hole in this radio 
mode is described by: 



Mi 



BH,R = K AGN 



/ M 



BH 



V10 8 M Q 



/hot 

~oT 



200 km s~ 



(7) 

where /hot is the fraction of the total halo mass in 
the form of hot gas and kagn is a free parameter 
controlling the efficiency of accretion. A value of 
kagn = 7.5 x 10~ 6 A1© yr -1 is found to reproduce 
the turnover at t he bright end of the galax y lumi- 
nosity function (|De Lucia fc Blaizodl2007t ). This 
simple phenomenological model may represent the 
accretion of cold gas clouds, or Bondi-Hoyle accre- 
tion fro m hot gas that fills the space between these 
clouds (jCroton et al.ll2006l ). 

The mechanical heating associated with radio 
mode accretion is assumed to impede, or even pre- 
vent, the cooling flow in central regions of the halo. 
More specifically, the radiated luminosity is taken 
to be 



LbK.R — £ r MBH,RC , 



(8) 



where c is the speed of light and e r describes how 
efficiently matter can be converted to energy near 
the event horizon. The standard value e r = 0.1 is 



adopted, which is an approximate value for radia- 
tively efficient accretion onto a non-rapidly spin- 
ning black hole (jShakura fc Svunvaevlli973f i. This 
injection of energy reduces the rate at which hot 
gas is able to cool onto the cold disk from M coo i 
to M coo i — 2LBH/^vir' wrtn the restriction that 
the cooling rate remains non-negative. The ef- 
fectiveness of radio mode AGN feedback in sup- 
pressing cooling flows is greatest at late times and 
for large black hole masses, precisely what is re- 
quired to reproduce the luminosities and colours 
of low-redshift, bright galaxies. 

2.3. Hydrodynamical simulations 

To explore the effect of energy feedback from 
galaxies on the properties of the ICM, we couple 
the L-Galaxies SAM to hydrodynamical simula- 
tions. The initial conditions for these simulations 
are the same as those used for the two dark matter 
simulations described above, but we add gas par- 
ticles with zero gravitational mass. The gas par- 
ticles then act purely as 'tracers' of the dark mat- 
ter. In this way we ensure that the dark matter 
distribution remains unaffected by the inclusion of 
baryons, so that the halo merger trees used to gen- 
erate the semi-analytic galaxy catalogues are the 
same. 

Recent work has shown that the dissipative 
nature of the baryon fluid ca n have an inf l uence 
on the structure of haloes dStanek et alj 120091: 
Romano-Dfaz et al.ll2009tlPedrosa et al.l2009r) an d 
the growth of merger trees (jSaro et alJ 12008). 
However, all semi-analytic modelling to date ig- 
nores such complications and uses merger trees 
based solely on the dark matter distribution. We 
are forced to follow this route so that we can use 
SAM input into our simulations. For the pur- 
poses of the results presented in this paper, the 
dark matter approximation is unlikely to make 
any significant difference to our conclusions. 

The number of gas particles that we add is 
iV gas = 100 3 for our small [L = 62.5/i^ 1 Mpc) 
simulation volume, and N gas — 200 3 for our larger 
(L = 125/i _1 Mpc) volume. We choose to in- 
clude gas at (approximately) the lower resolution 
of the Millennium Gas simulations H so that out 



2 The Millennium Gas simulations (Pearce et al., in prepa- 
ration) arc a suite of large hydrodynamical simulations, 
all having the same volume as the Millennium simulation 
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technique remains computationally feasible when 
applied to the full Millennium volume in future 
work. In Section [4] we conduct a resolution test to 
demonstrate that increasing the resolution of the 
gas component relative to the dark matter has a 
negligible effect on our results. 

The initial conditions are evolved from z\ = 127 
to z = with a modified version of GADGET- 
2, designed to allow for gas particles with zero 
gravitational mass. Whenever a SPH calcula- 
tion is to be done, we assume flb.o — 0.045 and 
assign the gas particles their corresponding true 
mass: m gas = 3.05 x 10 9 h~ 1 M@. This guarantees 
that gas properties such as density and entropy 
are computed correctly. In addition, gas parti- 
cles are also given their true mass for simulation 
data dumps, with the mass of the dark matter 
particles accordingly reduced to (1 — /b) m DM = 
7.05 x I0 8 h- 1 M Q , where f b = fi b ,o/fi ro ,o is the 
mean cosmic baryon fraction. 

In all of our hydrodynamical simulations we 
choose to neglect gas cooling processes since ra- 
diative cooling, star formation, black hole growth 
and associated feedback are incorporated in the 
SAM. Note, however, that gas particles are still 
converted to dissipationless 'star' particles as dic- 
tated by the SAM, following the procedure out- 
lined in Section [2.4. II below. 

We appreciate that the distribution and cool- 
ing of gas in haloes is treated in rather a simplistic 
manner in existing SAMs. However, we emphasise 
that, to the best of our knowledge, this is the first 
time the effect of energy feedback from galaxies 
on the ICM has been investigated using a SAM 
coupled to hydrodynamical simulations. As such, 
we believe it is sensible to begin with the simplest 
possible model where cooling is driven entirely by 
the SAM. In any case, the inclusion of gas cool- 
ing could only lead to a decrease in the entropy of 
intracluster gas, which would reinforce our conclu- 
sions with regard to the high degree of feedback 
that is required to explain the high entropy levels 
found in clusters. 



(L = 500/i — 1 Mpc) and utilising the same amplitudes and 
phases for the initial perturbations. The cosmological pa- 
rameters are also identical, except that the present baryon 
density parameter Cb,o = 0.045 to reflect the inclusion of 
gas. The simulations contain iVpjM = 1000 3 /2 dark matter 
and N gSbS = 1000 3 /2 gas particles, with respective masses 
1.42 X 10 w h- 1 M Q , and 3.12 X 10 9 h~ 1 M s . 



A more self-consistent approach would be to 
include radiative cooling in our hydrodynamical 
simulations and use the gas distribution to in- 
form the SAM. This extension of the semi-analytic 
technique would require the simulation and the 
SAM to be coupled in such a way that both can 
be undertaken simultaneously. Extensive testing 
would be necessary to ensure that such a model 
was as successful as current SAMs in reproduc- 
ing observed galaxy properties. Such a scheme is 
a long-term goal of our work but is beyond the 
scope of this paper. 

The semi-analytic galaxy catalogues that we 
have generated contain the positions and prop- 
erties of all model galaxies at 64 redshift values, 
corresponding to the output times of our simula- 
tions. We have further modified GADGET-2 so 
that, once an output redshift is reached, tempo- 
rary 'galaxy' particles are introduced at the appro- 
priate locations in the simulation volume. These 
galaxy particles have a set of associated properties, 
such as the change in stellar mass and energy re- 
leased by SNe/AGNs since the last output, which 
are calculated from the SAM galaxy catalogues 
prior to the simulation. We use this information 
to form stars and heat gas in the neighbourhood of 
each galaxy, as described in detail in the next sec- 
tion. Following star formation and the injection of 
energy, the galaxy particles are removed and the 
simulation progresses until the next output time, 
when the process is repeated. 

Note that the properties of the ICM should not 
be affected by the frequency with which energy is 
injected as the time interval between our chosen 
64 model outputs is always less than the galaxy 
halo dynamical time t^ yn — r v i T /v v i T = 0.1/H(z). 
We have verified that increasing the temporal res- 
olution has a negligible effect on our results. 

Cluster catalogues are constructed at z — for 
our simulatio ns using a procedu r e sim ilar to that 
employed by iMuanwong et al. ( 2002 ). Briefly, 
groups of dark matter particles are identified with 
the FOF algorithm, setting the linking length to 
be 10% of the mean inter-particle spacing. A 
sphere is grown about the most gravitationally- 
bound dark matter particle of each group until 
radii are found that enclose mean overdensities 
of A = 94, A = 200, A = 500, A = 1000 and 
A = 2500, relative to the critical density p CI} Q. 
Any clusters which overlap with a more massive 



S 



system within these radii are discarded from the 
catalogues. 

To check that dark matter structures are in- 
deed undisturbed by the massless gas particles, 
we have compared the positions of cluster centres 
in catalogues generated from our dark matter-only 
simulation in the L = 62.5/i _1 Mpc box and a non- 
radiative hydrodynamical simulation in the same 
volume. At each overdensity A, we find that the 
distance between corresponding cluster centres in 
the two catalogues is less than the softening length 
for over 99% of our identified objects. We have 
also verified that this is the case at high redshift 
z » 3. 

2.4. Implementing star formation and 
feedback from galaxies 

2.4-1- Star formation 

Each model galaxy in the semi-analytic cata- 
logues has an associated stellar mass M* . To com- 
pute the mass locked-up in stars formed over a red- 
shift interval Az — z n+ \ — z n , we take the stellar 
mass of the galaxy at z n and subtract the sum of 
the stellar masses of its progenitors at the previous 
output Z n+ l- 

AM, = M*0„) - M *i z n+i)- (9) 

prog. 

Once an output redshift is reached in a simu- 
lation, we convert the AJV 8 t ar = AM,/ro gas gas 
particles nearest to each model galaxy into col- 
lisionless star particles, reflecting the increase in 
stellar mass AM, of the galaxy since the last out- 
put. To ensure that AiV s t a r is an integer, we draw 
a random number r uniformly from the unit in- 
terval and compare it with the fractional part of 
AJVgtar: if r is less (greater) than the fractional 
part of AJVgtari we round AjV s t ar up (down) to the 
nearest integer. Note that the star particles also 
have zero gravitational mass, so they do not influ- 
ence the dark matter distribution. 

2-4-2. Type II supernova feedback 

The energy released by Type II SNe between 
two successive outputs depends on the amount of 
mass that went into new stars during this time 
period. For each model galaxy, this is simply 
Am, = AM,/(1 — R), where AM, is obtained 



from equation ((9|). Inserting Am, into equation 
([5|) gives the mass of ejected gas, corresponding to 
an energy input into the ICM of 

A£ CJCCtcd = iAm ej eoted<4 (10) 

over the interval Az. The solid line in Fig. [T] 
shows the L-Galaxies prediction for the cumulative 
total energy transferred to the ICM by Type II 
SNe as a function of redshift. For comparison, the 
maximum possible total energy available from SNe 
is also shown as the dotted line. This is computed 
at each redshift by cumulatively summing A_EgN 
for all galaxies, assuming an efficiency ehaio = 1< 

At each output, we distribute the available en- 
ergy Adjected amongst the neighbouring gas par- 
ticles of a galaxy using a simple heating model. 
The basis of our model is that all gas within a 
distance r V i r of a galaxy can be heated by feed- 
back processes in the time At between two out- 
puts. We choose to inject energy in a distributed, 
rather than local, manner since, over a time At, 
heated gas will flow outwards, mixing with in- 
falling cooler gas at larger radii. We have also ex- 
perimented with alternative heating models where 
each galaxy heats a fixed number of neighbouring 
gas particles and found that our results are not sig- 
nificantly affected by the choice of heating model. 
In a forthcoming paper we intend to investigate 
different ways of injecting energy into the ICM in 
more detail. 

We implement our heating model in GADGET- 
2 as follows. For each galaxy, we find the number 
Angb °f S as particles within a sphere of radius r v ; r 
centred on the galaxy. If no neighbours are found, 
the search radius is increased until one gas particle 
is found. This is typically only necessary for low 
mass haloes. The energy released by SNe since the 
last output is then used to raise the entropy of the 
neighbouring gas particles by a fixed amount 

AA l = (7 -Rejected (n) 

™ g as E,=f [max (/ bPvir , Pj )\i-i 

where A = /cT/(/im p p 7_1 ) is the definition of en- 
tropy employed by GADGET-2, m p is the mass 
of a proton and /i w 0.6 is the mean molec- 
ular weight for a fully-ionised gas of primordial 
composition. By giving each particle a fixed en- 
tropy, rather than energy, boost, we ensure that 
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Fig. 1. — Cumulative total amount of energy 
transferred to the ICM by different feedback 
sources as a function of redshift in a box of side 
length L = 62.5ft." 1 Mpc. The dotted line shows 
the total energy released by supernova explosions, 
while the solid line is the fraction of this energy 
available to heat intracluster gas. The contribu- 
tions from quasar mode AGN feedback when the 
thermal coupling efficiency is ef = 1, ef = 0.4 and 
£f = 0.05 are given by the dashed, dot-dashed 
and long dashed lines, respectively. The triple- 
dot-dashed line shows the mechanical heating as- 
sociated with radio mode accretion. However, in 
the Munich L-Galaxies semi-analytic model, this 
energy only reduces the rate at which gas can cool 
out of the hot halo, rather than heating the ICM. 



denser particles close to the galaxy are heated to a 
higher temperature than more distant, lower den- 
sity particles. The product of the cosmic baryon 
fraction /b and the virial density p V ir gives the 
mean overdensity of baryons within the virial ra- 
dius. If no neighbours are found within a dis- 
tance r vir of a galaxy and the search radius has to 
be increased until one particle is found, the den- 
sity of this particle may be less than /b/o v ir- By 
using [max (/bPvir, Pj)] 7 " 1 , rather than Pj~ , in 
the sum in equation (fTTj) . we are assuming that 
the amount of energy used to heat such particles 
is A£ C j CCt od(Ai//bPvir) 7 ~ 1 < A.E e je C ted; the rest 
of the energy is taken to be used up as the gas 
does work expanding adiabatically to a density 

Pi < /bPvir- 

2.4-3. Feedback from Active Galactic Nuclei 

A similar approach is used to calculate the en- 
ergy released by AGN activity over a redshift in- 
terval Az. For each galaxy, the total change in 
mass of the central black hole AMbh is given by 
equation ([9]), but with M* replaced by Mbh- The 
fraction of this mass accreted via the radio mode 
is approximately 

AM BH ,r - M BH<R At, (12) 

where Mbh,r is obtained by evaluating equation 
([7|) at the current output z n . It follows that 
the mass change due to merger-driven accretion 
is then 

AM BH ,q = AM BH - AAf BH ,R. (13) 
We ensure that the radiated luminosity 

Lbh.q = e r M B H,QC 2 (14) 
does not exceed the Eddington luminosity 



e.g. Begelman fc Meierll982 ) when averaged over 



the time At between two successive outputs. We 
then assume that some fraction of the energy ra- 
diated during quasar mode accretion is coupled 
thermally to the ICM: 
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(16) 



where the coupling efficiency ef is a free parame- 
ter; we discuss suitable values for this parameter 
in the following section. The dashed, dot-dashed 
and long dashed lines in Fig. [T]show the cumula- 
tive total energy transferred to the ICM by quasar 
mode AGN feedback for the cases ef = 1, ef = 0.4 
and ef = 0.05, respectively. We inject the energy 
ALbh.q released by quasar mode accretion into 
gas particles surrounding model galaxies using the 
same heating model as for supernova feedback. 
The mechanical heating 



A£bh,r = e r AA/ B H,RC 2 



(17) 



associated with quiescent accretion simply reduces 
the rate at which gas in the hot halo cools onto 
the disk. In other words, radio mode feedback 
does not explicitly heat intracluster gas within 
the framework of L-Galaxies and is thus irrelevant 
for our hybrid approach. For completeness, we 
show the cumulative total energy liberated by ra- 
dio mode accretion as the triple-dot-dashed line in 

Fig. m 

3. RESULTS AND DISCUSSION 
3.1. Model discrimination 

One of the most fundamental properties of 
galaxy groups and clusters is the X-ray luminosity- 
temperature relation. Only recently have self- 
consistent hydrodynamical simulations been able 
to successfully reproduce the slope and normali- 
sation of this relation over a wide range of mass 
scales. In this section we explore the individual 
effects of stellar and AGN feedback on the Lx-T 
relation with a set of test simulations performed 
in the L = 62.5/i~ 1 Mpc box using our hybrid 
technique. 

3.1.1. Star formation and supernova feedback 

As a starting point, we include only star for- 
mation and associated supernova feedback as pre- 
dicted by L-Galaxies. Fig. [2] shows our simu- 
lated groups and clusters on the Lx-T plane at 
z = 0, shaded by their hot gas fraction. Bolo- 
metric X-ray luminosities and emission-weighted 
temperatures were computed using the procedure 



outlined by iMuanwong et al ( 2002 ). Since radia- 
tive cooling is neglected in our simulations, we do 
not anticipate excess emission in cluster cores due 
to strong cooling flows and it is therefore unnec- 
essary to remove emission from central regions. 
All cluster properties are calculated within rsoo 
since A = 500 is typically the smallest density 
contrast accessible to observations. Furthermore, 
we only consider galaxy groups and clusters with 
a mass Af 50 o > 1.3 x lO 13 ft. _1 Af0, correspond- 
ing to a total of about 15,000 particles within 
rsoo. The most massive object formed has a mass 
M 500 w 7 x l0 13 hr 1 M Q . 

For comparative purposes, data from a num- 
ber of X-ray observational studies is shown 



(lMarkevitchll998tlArnaud fc Eyrardll999tlWu et al 
I l999t iHelsdon fc Ponmaij l200d: IXue fc Wull200oT~ 



Horner 2001 uMulchaev et a l. 2003; Osmond fc Ponman 
20041 ). The uppermost solid line in Fig. [5] is the 
best-fit Lx~T relation obtained from a simulation 
that includes gravitational heating only. We find 
that Lx oc j ,1 - 91±0 - 37 in this case, close to the 
Lx oc T 2 scaling predicted by theoretical argu- 
ments. The lower solid line shows the best-fit 
Lx~T relation for a simulation with uniform pre- 
heating at high-redshift. The preheating model 
adopted imposes an entropy floor <Sp re heat — 200 
keV cm 2 at z = 4. Preheating leads to a much 
steeper relation: Lx oc T 3 - 87±0 - 87 , consistent with 
observational data. 

Observe that nearly all of our data points lie 
just below the best-fit Lx-T relation for the grav- 
itational heating simulation, with only a slight 
hint of steepening. The hot gas fractions of our 
groups and clusters are at least 65% of the cos- 
mic baryon fraction, with some having a hot gas 
fraction very close to /b. Observational data sug- 
gests that only massive clusters (T > 5 keV) have 
such large gas fractions, with groups and poor 
cluste rs typically having much smaller gas frac- 



tions (ISanderson et al.ll2003tlVikhlinin et al.| [2006: 



Sun et al J 120091) . These results indicate that the 
entropy of the ICM has been raised by stellar 
feedback, but nowhere near enough to explain ob- 
served X-ray luminosities. 

By contrast, hydrodynamical simulations with 
radiative cooling, star formation and supernova 
feedback tend to produce a closer match to the 



observed Lx-T relatio n (e.g. iBorgani et al.l 12004 ; 
Puchwein et al. 2008). There are two reasons 
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Fig. 2. — Bolometric X-ray luminosity as a func- 
tion of emission-weighted temperature for groups 
and clusters formed in a hybrid simulation with 
L-Galaxies star formation and Type II supernova 
feedback. X-ray properties are computed within 
r5on and each object is shaded according to its 
hot gas fraction. For comparison, the upper and 
lower solid lines are the best fit relations from non- 
radiative and preheating simulations, respectively. 
Several datasets obtained from X-ray observations 
of groups and clusters are also shown. 



for this. Firstly, the amount of baryons that 
cool and form stars is typically much greater in 
hydrodynamical simulations, with 30 — 50% of 
the baryons within the virial radius locked-up in 
stars (e.g. IBalogh et all |200ll iDave et al.l [2002; 



Tornatore et alj l2003h . In our hybrid simulation 



the average stellar fraction within the virial radius 
is approxima tely 9%, in agreement with observa- 



tional data (IBalogh et all l200ll iLin et all [2003 



Balog h et al 



2008). This is, of course, to be ex- 



pected since star formation in our simulation is 
driven by a SAM which has been tuned to re- 
produce the cosmic star formation history. Sec- 
ondly, popular stellar feedback schemes typically 
assume that each supernova event transfers more 
energy to the surrounding gas. For e xamp le, in 
the model s of ISpringel k, Hernquistl (|2003l ) and 
Kay et al. ( 2004f i!. the amount of energy returned 



to the ICM per solar mass of stars formed is at 
least 4 x 10 48 erg Mg 1 . Within the framework of 
L-Galaxies, it follows from equation ^ that this 
is only possible if no energy is used to reheat cold 
disk gas (AiShot = 0) and the efficiency parameter 
Ehalo = 1- 

To test that our method yields comparable re- 
sults to existing hydrodynamical simulations when 
the star formation efficiency is increased, we have 
performed a simulation where the mass of newly 
formed stars Am* is simply multiplied by 5. We 
also assume the energy imparted to the ICM by 
SNe is larger: A_E C j OCtc d = Ai?gN, with an effi- 
ciency ehaio = 1- Note that setting Ai?hot = 
means all star-forming galaxies can eject gas from 
their halo, rather than just those with a virial ve- 
locity w V ir ^ 200 km s _1 . Fig. [3] shows the re- 
sulting Lx~T relation. It is evident that cluster 
X-ray luminosities and hot gas fractions are sig- 
nificantly reduced relative to those predicted by 
the standard L-Galaxies stellar feedback model. 
Our results are in broad agreement with those ob- 
tained from simulations incorporati ng star forma- 



tion and supernova feedbac k (e.g. iBorgani et al 



2004 IPuchwein et al.l l2008h . However, the nor- 



malisation of the Lx-T relation is still too high and 
the slope is too shallow relative to the observed re- 
lation. The stellar fraction within the virial radius 
has now increased to 39% on average, conflicting 
with observations but a gain similar to results from 
direct simulations (e.g. Puchwein et al. 20081 ). 
The heating of intracluster gas by stellar feed- 
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back from model galaxies is clearly insufficient 
to reproduce the Lx-T relation, particularly on 
group scales. This remains the case even when 
the star formation efficiency is unrealistically high. 
Even if we were to allow for the energy released 
by Type la SNe, this is unlikely to bring the Lx-T 
relation shown in Fig. [2] in line with observations 
since the energetics of Type la and Type II SNe are 
thought to be roughly similar. However, Type la 
SNe are crucially important for c hemically enrich 



ing t he ICM at low-redshift (e.g. iMannucci et al 
2005); an issue we shall address elsewhere. We 



now investigate whether additional energy input 
from AGNs can resolve this problem, starting with 
feedback from quasar-induced outflows. 

3.1.2. Including feedback from Active Galactic 
Nuclei 

Cosmological simulations including prescrip- 
tions for black hole growth and associated quasar 
activit y have been performed bv lDi Matteo et aT 
(200&j). They assumed 5% of the energy radiated 



during quasar mode accretion is coupled thermally 
and isotropically to the surrounding medium. 
This choice ensures that the normalisation of the 
local black h ole-bulge mass re l ation agrees with 



observations ()Di Matteo et al.l 120051 ). In their 
model, the radiated luminosity is related to the 
total black hole accretion rate Mbhj whereas we 
assume it is governed by the merger-driven accre- 
tion rate Mbh,q- However, the growth of black 
holes in L-Galaxies is dominated by quasar mode 
accretion, so that AA/bh ~ AA/bh, q anyway. 
Based on this argument, it seems sensible to be- 
gin by setting et = 0.05 in equation (fTC|) . 

Fig. [4] shows the Lx-T relation obtained from 
a hybrid simulation including star formation, su- 
pernova and quasar mode AGN feedback with a 
coupling efficiency ef = 0.05. Again, cluster X-ray 
luminosities and hot gas fractions a re con siderably 



over-estimated. |Pi Matteo et al. (2008) did not 



use their model to investigate the Lx-T relation 
of groups and clusters. However, since the global 
black hole accretion rate histories in L-Galaxies 



and the simulations of iDi Matteo et al.l (|2008f ) are 
similar, we would expect them to obtain a similar 
result. 

In order to recover the observed steepening of 
the Lx~T relation on group scales, we find that 
the quasar mode coupling efficiency has to be in- 
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Fig. 3. — Same as Fig. [2 but with an enhanced 
stellar feedback scheme in which the L-Galaxies 
star formation efficiency is 5 times greater and the 
amount of energy transferred to the ICM by Type 
II SNe has been set to its maximum possible value. 
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Fig. 4. — Same as Fig. [2j but for a hybrid simula- 
tion with L-Galaxies stellar feedback and a simple 
model for quasar mode AGN feedback. The energy 
released by merger-driven accretion onto the cen- 
tral black hole is assumed to be thermally coupled 
to the intracluster gas with an efficiency £f = 0.05. 
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creased to ef = 0.4. The Lx-T relation in this 
case is shown in Fig. [5] Our simulated clusters 
now provide a much better match to the observa- 
tional data. In addition, the gas fractions of our 
clusters are in broad agreement with ob servations 



(jSanderson et alJl2003t[Sun et al.ll2009l) 



Artificially boosting the feedback from merger- 
driven accretion can be viewed as a crude attempt 
at including radio mode feedback since, at each 
redshift, over 75% of model galaxies with a central 
black hole accrete via both modes. Can we recover 
the correct Lx-T relation by explicitly including 
the radio mode contribution, rather than simply 
increasing ef? Recall that the L-Galaxies imple- 
mentation of radio mode feedback only reduces the 
rate at which gas can cool out of the quasi-static 
hot halo, rather than directly heating the ICM. 
Even if we temporarily ignore this feature of the 
model and inject all of the energy AEbh.r (see 
equation I17[) released by quiescent accretion into 
the gas component of a simulation incorporating 
stellar and quasar mode feedback with ef = 0.05, 
we find that it has little effect on the Lx-T rela- 
tion shown in Fig. [4] This is because the mechan- 
ical heating associated with radio mode accretion 
is negligible at high redshift, only becoming com- 
parable to the energ y released by SNe at z = 0; 
see Fig. [T] However, iPuchwein et alJ (|20081 ) have 
successfully reproduced the observed mean Lx-T 
relation by incorporating a model for AGN-driven 
bubble h eating into the quasar mode feedback 
scheme of l5i Matteo et~aH (j2008l) . This highlights 
the need for an alternative model of AGN feedback 
in which the radio mode does more than just offset 
cooling from the hot halo. 

In recent work, iBower et al. (2008) have used 
the Durham SAM GALFORM to investigate the 
properties of the ICM, particularly the Lx-T re- 
lation. They demonstr ated that the vers ion of 
the model developed by IBower et al. (200G) over- 
predicts X-ray luminosities on group scales, lead- 
ing to a shallow Lx-T relation similar to the one 
we obtained from our simulation with L-Galaxies 
stellar feedback (Fig. [2]). In this model, the energy 
released by quiescent radio mode accretion simply 
prevents any significant amount of gas cooling in 
massive haloes, as in L-Galaxies. To try to explain 
the X-ray properties of galaxy groups and clus- 
ters whilst simultaneously ac counting for the ob - 
served properties of galaxies. IBower et al. (2008) 



proposed a modification of GALFORM which al- 
lows for heat input into the ICM from radio mode 
feedback. This additional heating acts to expel 
gas from the X-ray emitting central regions of 
haloes, reducing the gas density and thus luminos- 
ity. Lower-mass systems arc affected more than 
massive ones because their cooling time is shorter, 
meaning they initially supply more material to the 
central black hole, resulting in a larger amount of 
feedback per unit mass of gas. The Lx-T relation 
then becomes steeper as desired on group scales. 
With this modification, the GALFORM model re- 
produces the observed mean Lx-T relation and 
the substantial scatter about this relation at low 
temperat ures (T < 3 keV) , as well as halo gas 



fractions ( Bower et al. 200 



The basis of the lBower et all (|2008f ) AGN feed- 
back model is to compute the heating power Lheat 
as a function of the cooling rate. More specifically, 



-^hcat = ^SMBHCr-McoolC 2 , 

subject to the constraint 

Lhcat < CSMBH^Edd- 



(18) 



(19) 



Here M coo i is the cooling rate of the halo in the 
absence of radio mode feedback, related to the 
black hole accretion rate via Mbh = ??smbh^coo1- 
The parameter 7/smbh controls the efficiency 
with which cooling material can be accreted by 
the black hole and is set to jjsmbh = 0.01 by 



Bower et al 



(|2008l) . The limiting criterion (TT 
relates to the structure of the accretion disk itself. 
Effective radio mode feedback requires efficient 
jet production, typically thought to be associ- 
ated with geomet rically thick, advection domi- 
nated disks (e.g. iRees et all 119821 iMeierl 12001 
Churazov et al.ll2005l ). If the accretion rate is too 
high, models suggest that the vertical height of the 
disk will collapse, leading to a drop in jet produc- 
tion efficiency. Much more of the energy released 
by accretion is then radiated away and is not 
available for radio mode feedback. Based on the- 



oretical work (e.g. lEsin et al.lll997l ). IBower et al 



([20081 ) assume this change in structure of the ac- 
cretion disk occurs once the accretion rate reaches 
M B h = esMBH^Edd, where csmbh is referred to 
as the disk structure parameter and is related to 
the disk viscosity parameter a by esMBH °c o 2 - 
The disk structure parameter is assigned a value 
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esMBH = 0.02, in accordan ce with plausible accre- 
tion disk viscosities (e.g. iMcK innev fc Gammiel 
l2004t iHirose et alJ l2004 lHawley fc Krohkll2006h . 
Note that th e AGN feedback model suggested by 
Siiacki et all (|2007t) also assumes that radio mode 



feedback is only effective for accretion rates below 
some fraction of the Eddington rate. 

The heating energy available from the AGN is 
compared with the energy lost radiatively by gas 
cooling from the halo. If the feedback energy is 
greater than the radiated energy, the excess energy 
is used to eject gas from the halo. The energy 
transferred to the ICM over a time period At is 
then 



A-E e jectcd — AEfrcat — AE coo \, 

where AiJhcat = L heat At and 



AE. 



1 



cool 



AM cool 4. 



(20) 



(21) 



We have implemented the Bower et al.l ( 20081 ) 
prescription for AGN feedback in our hybrid ap- 
proach using the following procedure. For each 
galaxy in the catalogue, we take the total change 
in black hole mass between two successive outputs, 
AM B h, and compute AAf cool = AA/ B h/77smbh- 
Inserting this into equations (fT5| and ([2"!]) gives 
A-Ehcat and AE coo \, respectively. The heat input 
into the ICM from each model galaxy then fol- 
lows from equation (|20p. This energy is injected 
into the gas component of our simulation using the 
same heating model as described previously. 

To be fully consistent, we should take into ac- 
count the effect of AGN heat ing on hot halo ga s 
within L-Galaxies itself. As in lBower et al" (2008), 
this would almost certainly lead to changes in 
some of the other model parameters. However, if 
we were to adjust the appropriate parameters so 
as to regain a realistic galaxy population, then we 
would require galaxy properties (particularly their 
black hole masses) to be similar to those in our 
existing galaxy catalogue. The amount of heating 
energy available from AGNs would then also be 
(approximately) the same since it is driven purely 



by the black hole growth rate in the iBower et al 
(2008) model. 



Fig. [6] shows the Lx-T relation obtained from 
a simulation with L-Galaxies stellar feedback and 



the AGN feedback model of IBower et al.l (|2008fi . 



Our results clearly provide an excellent match to 
the observational data. In addition, the hot gas 
and stellar fractions of our groups and clusters 
agree with observed values. This success indicates 
that the L-Galaxies implementation of AGN feed- 
back must indeed be revised if the model is to ex- 
plain the observed properties of both the galaxy 
distribution and the ICM. 

3.2. Best-fit model 

Following the analysis presented in the preced- 
ing section, we have established a model that is 
capable of reproducing the observed Lx-T relation 
of groups and poor clusters. This best-fit model 
comprises of the stellar feedback scheme employed 
in L-Galaxies, combined w i th the AGN feedback 
prescription of Bower et al. ( 20081 ). To investigate 
whether this model can also correctly explain the 
properties of richer clusters, we have performed 
a hybrid simulation in a larger L — 125/i _1 Mpc 
volume. 

3.2.1. The X-ray luminosity-temperature rela- 
tion 

The z = Lx~T relation obtained from our 
simulation is shown in Fig. [7j Groups and clus- 
ters are shaded by their hot gas fraction and we 
only plot clusters with a mass M500 > 1.3 x 
10 13 /i _1 M©. The mass of the largest cluster 
formed is now M 500 ~ 2.3 x 10 u ft- 1 M s . For 
comparison, we show the best-fit Lx-T relations 
from a non-radiative simulation (upper solid line) 
and a simulation with the same preheating model 
as described previously (lower solid line). These 
relations are of the form Lx oc T 1 - 88 ^ - 16 and 
Lx oc T 3 62±0 - 33 , respectively. As before, our hy- 
brid simulation yields an Lx-T relation with a 
slope and normalisation that is generally consis- 
tent with observations at all mass scales. A few 
objects appear to have a higher temperature than 
expected, given their X-ray luminosity. This is 
particularly true for the T ss 3.3 keV cluster with 
a luminosity Lx ~ 5.4 x 10 39 erg s _1 . The reason 
for this is that the gas in this object has recently 
been raised to a high temperature by a large en- 
ergy injection from the central AGN, causing it 
to flow outwards. As this gas is replaced by in- 
falling cooler gas, the system will stabilise and 
shift back towards the main relation. Our results 
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Fig. 5. — Same as Fig. IH except the thermal 
coupling efficiency for quasar mode AGN feedback 
has been increased to ef = 0.4. 
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Fig. 6. — Same as Fig. [2J but for a hybrid 
simulation with L-Galaxies stellar feedback and 



the Bower et all (j2008f ) AGN feedback prescrip- 
tion employed in the latest version of GALFORM. 



exhibit a variation in scatter along the Lx-T re- 
lation that is similar to observational data. In 
particular, for temperatures T < 3 keV, our data 
points fan out to populate a triangular region of 
the Lx~T plane in the same way. This scatter 
about the low-temperature end of the Lx-T rela- 
tion is attributable to the variety of merger histo- 
ries of groups. Note that there seems to be more 
scatter towards the lower-luminosity edge of the 
observed relation than the upper-luminosity edge. 
This is because we cannot produce systems with 
a highly X-ray luminous CC in our simulations 
since we have chosen to neglect cooling processes 
(see Section l3.2.3j) . 

3.2.2. Halo gas fractions 

It is evident from Fig. [7]that, as we move along 
the Lx-T relation from group to cluster scales, the 
hot gas fraction of our simulated objects increases, 
reaching approximately 70% of the mean cosmic 
baryon fraction in massive clusters. To illustrate 
this behaviour more clearly, we explicitly plot the 
hot gas fraction / gas as a function of emission- 
weighted temperature (both within r5oo) in Fig. 
IU For comparative purposes, we show constraints 
on hal o gas fractions obtained from X-ray ob serva- 
tions (|Vikhlinin et alJl200& ISun et alJl2009h . Also 
plotted are gas fractions within rsoo computed 
from the gas den sity and tempera t ure p rofile pa- 



rameters given in ISanderson et al. I (|2003l) 



The gas fractions of our groups and clusters are 
in broad agreement with the observational data, 
exhibiting a comparable amount of scatter. In 
both our results and the data we see a rapid de- 
cline in gas fraction at lower temperatures. This is 
because AGN feedback is effective at driving gas 
from the central regions of low-mass systems. On 
the other hand, the potential wells of massive clus- 
ters are too deep for AGN heating to efficiently re- 
move gas from them, so they retain a much larger 
fraction of their hot gas. Finally, we note that the 
T w 3.3 keV cluster with the smallest gas frac- 
tion (0.3%) corresponds to the system in Fig. [7] 
with Lx ~ 5.4 x 10 39 erg s _1 discussed before. 
This supports our argument that this object has 
recently experienced an intense burst of AGN ac- 
tivity which has driven a substantial amount of 
gas beyond r 500 . 
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3.2.3. Entropy and temperature profiles 
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Fig. 7. — Bolomctric X-ray luminosity as a func- 
tion of emission- weighted temperature for a hybrid 
simulation with our best-fit feedback model. This 
model consists of s tellar feedback from L-Galaxies 
combined with the lBower et~ail (|2008f) AGN feed- 
back prescription. See the caption of Fig. [2] for a 
description of the shading and best-fit lines. 
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Fig. 8. — Halo hot gas fraction within rsoo as 
a function of emission-weighted temperature for 
a hybrid simulation with our best-fit feedback 
model. Simulated objects are shown by the open 
circles. A variety of data from X-ray observations 
of groups and clusters is shown for comparison. 
The horizontal dashed line is the mean cosmic 
baryon fraction, /b, in the cosmological model we 
have adopted. 



Further information about non-gravitational 
heating processes operating in clusters can be 
gleaned by inspecting radial entropy and temper- 
ature profiles of the ICM. In the upper (lower) 
panel of Fig. [51 the thin solid lines are en- 
tropy (emission-weighted temperature) profiles 
for a sample of five of our most massive clus- 
ters, with masses in the range 1.1 x 10 14 /i -1 M© < 
M 500 < 1.7 x IO^/i^Mq at z = 0. Thick solid 
lines highlight the profiles of the most massive 
cluster in this sample (the second most massive 
object in the simulation volume). For reference, 
the thick dotted and dashed lines are the pro- 
files of the corresponding object formed in the 
non-radiative and preheating simulations, respec- 
tively. Profiles are only plotted for radii greater 
than the gravitational softening length. In addi- 
tion, observed profiles for four NCC cluste r s from 
the sample presented by ISanderson et al. (2009) 
are shown as thin dot-dashed lines. These ob- 
jects are of similar mass to our simulated clusters: 
1.1 x IO^/i^Mq < Af 50 o < 1-6 x 10 14 ft- 1 M Q , 
with the profiles of the most massive object high- 
lighted by thick dot-dashed lines. We only com- 
pare our profiles with those of NCC clusters since 
our simulations are non-radiative and thus sys- 
tems with a CC do not form. 

Theoretical models of shock heating during 
spherical collapse predict tha t entropy scal es 



with radius as S oc r (e.g. iTozzi fc Norman 



2001). This scaling behaviour is indeed observed 
in cluster outskirts, but entropy profiles are typ- 
ically seen to b ecome flatter in central regions 



r < .2r 2 on (e.g. Ponman et al.l 2003 ; Pratt et al 
2006). However, the precise radius at which this 
flattening occurs varies considerably, depending 
on such factors as the temperature (mass) of the 
system and whether it has a CC or a NCC. In 
particular, hotter, more massiv e objects have a 
highe r mean core entropy (e.g. ICavagnolo et al. 
2009), and the profiles of NCC clusters flatten 



off at signific antly larger radii than those of CC 
clusters (e.g. Sanderson et al. 20091 ). 

The power-law S oc r 11 is illustrated by the 
thin dashed line in the upper panel of Fig. \§\ (the 
normalisation is arbitrary). It is evident that the 
entropy profiles of our clusters depart from this 
scaling at somewhat larger radii than observed: 
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Fig. 9. — Radial profiles of entropy (upper panel) 
and emission-weighted temperature (lower panel) 
obtained from a hybrid simulation with our best- 
fit feedback model. The solid lines are profiles 
for a sample of five clusters with a mass 1.1 x 
IO^/i^Mq < M 500 < 1.7 X lO 14 ft- 1 M , with the 
profiles of the most massive cluster in this sample 
highlighted by thick solid lines. For comparative 
purposes, the thick dotted and dashed lines are 
the profiles of the corresponding object formed in 
a non-radiative and a preheating simulation, re- 
spectively. In addition, observed profiles for four 
no n-cool core clusters from the sample presented 
by ISanderson et all ( 20091 ) are shown as thin dot- 
dashed lines. These objects are of similar mass 
to our simulated clusters: 1.1 x lO^h^M™ < 



M 500 < 1.6 x lO u h -1 M , and the thick dot- 
dashed lines highlight the profiles of the most mas- 
sive object. Theoretical arguments predict that 
entropy scales as S oc r outside of cluster cores, 
shown by the thin dashed line in the upper panel. 



r ~ r5oo , flattening as we move in towards the core 
and over-estimating the central entropy. Likewise, 
the temperature profiles shown in the lower panel 
of Fig. [5]reveal that the core temperature is higher 
than would be expected for objects of this mass. 
This is similar to the behaviour predicted by the 
preheating model, although our feedback model 
yields an entropy profile that provides a slightly 
better match to the observational data. 

Based on the preceding discussion it is natural 
to ask how we can recover the observed Lx-T re- 
lation and halo gas fractions within r^o if the gas 
entropy is over-estimated in clust er cores. The an- 



swer lies in the fact that the full Sanderson et al 



(2009) sample contains clusters with a variety of 
entropy profiles. If we consider the emission per 
unit logarithmic radius, dLx/dlog 10 r, the most 
centrally-concentrated CC clusters have their peak 
emission within 0.2r 5 oo, whereas in some NCC 
clusters the peak is beyond 0.7r-5oo- Our most- 
massive simulated clusters resemble these latter 
objects, and hence lie at the lower-luminosity edge 
of the Lx-T relation (for low-mass systems, the ob- 
servations do not extend far enough in radius to 
make a meaningful comparison). 

We note that profiles can only be reliably mea- 
sured for groups and poor clusters that are very 
X-ray bright. However, optically-selected samples 
of groups have revealed systems with little or no 
detectable X-ray emission that may be the group 



Rasmussen et al 



analo gue of NCC clusters (e.g. 
20061 ). It may not be possible to extract profiles 



for such objects, but since they are X-ray under- 
luminous we would expect them to have a large 
entropy core, possibly similar to that seen in the 
profiles of our simulated groups. 

However, in order to reproduce the observed 
profiles of the majority of X-ray bright groups and 
poor clusters, the entropy and temperature of the 
gas in core regions must be lowered. This could 
potentially be achieved by incorporating cooling 
processes in our simulations. The inclusion of ra- 
diative cooling would only strengthen our conclu- 
sion that large amounts of energy must be injected 
into the ICM by AGNs to recover the observed 
Lx-T relation and halo gas fractions. 
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4. RESOLUTION TEST 

In all of the hydrodynamical simulations pre- 
sented in this paper we have chosen a lower res- 
olution for the gas component than for the dark 
matter. To check that this approach yields robust 
results, we have performed a suite of six hybrid 
simulations in the L = 62.5h~ 1 Mpc box with a 
fixed number TVdm = 270 3 of dark matter parti- 
cles, but different numbers of gas particles. Two of 
these simulations have N gas < 100 3 , while for the 
other three N gm > 100 3 (recall that all our simula- 
tions in this volume have N gaa = 100 3 ). All other 
simulation parameters are unchanged between the 
runs. The feedback scheme we adopt is our best- 
fit model since this is the only one that we have 
used in the larger L — 125/i~ 1 Mpc box and it is 
this model we intend to apply in the future. 

To demonstrate that N gaa = 100 3 is a sufficient 
number of gas particles for a converged estimate of 
the Lx-T relation, in Fig. [T0]we compare the rela- 
tions obtained from three simulations with N gas = 
50 3 , iV gas = 100 3 and A gas = 200 3 , respectively 
The corresponding (true) masses of the gas parti- 
cles are then m gas = 2.44 x 10 10 /i _1 Af Q , m gas = 
3.05 x lO^M©, and m gas = 3.81 x lO^M©. 
X-ray luminosities and emission-weighted temper- 
atures are computed within rsoo and only objects 
with a mass M 50 o > 1.3 x 10 13 /i -1 M Q are plotted. 
Focusing on the two most massive clusters, we see 
that they have very similar temperatures and lu- 
minosities in the intermediate and high-resolution 
runs. We find that this remains the case if we go 
to even higher resolution. However, they shift sig- 
nificantly in the low-resolution simulation. This 
signifies numerical convergence at the intermedi- 
ate resolution for systems with T > 2 keV, which 
is important since we plan to use our technique 
to preferentially study rich clusters in subsequent 
work. 

As we move to lower temperatures, the amount 
of scatter between the three different resolution 
simulations increases. In the low resolution run, 
the data points are less tightly grouped and seem 
to lie on a mean relation that has a shallower slope 
than in the other two simulations. In addition, 
three objects have anomalously high temperatures 
given their X-ray luminosity. This is because there 
is a limited number of gas particles available for 
heating in these low-mass systems at this resolu- 



tion and a large amount of energy has recently 
been injected into a small number of gas particles. 
On the other hand, the groups and clusters formed 
in the intermediate and high-resolution runs lie on 
a similar mean relation, although the data points 
do not line-up perfectly. This indicates that we 
have not attained numerical convergence at the 
lowest resolution. Note that the Lx-T relation in 
the intermediate and high-resolution runs appears 
to flatten slightly at the low temperature end in 
the same way as observational data. 

It is not surprising that we do not achieve exact 
convergence for two reasons. Firstly, in situations 
where only one gas particle is heated, it is evi- 
dent from equation (fTTj) that the entropy boost 
given to the particle is inversely proportional to 
its (true) mass: AAi oc l/m gas , if its density is 
less than /bPvir- Consequently as the resolution 
is increased, the particle will receive a larger en- 
tropy injection. Secondly, there is a stochastic ele- 
ment to our star formation scheme. By performing 
another intermediate-resolution simulation with a 
different random seed, we found that this random- 
ness induces considerable scatter in the luminosi- 
ties and temperatures of the lowest-mass systems, 
but has little effect in more massive objects. Based 
on these considerations, we suggest that the X- 
ray properties computed for groups with T < 1 
keV are unreliable since these systems contain rel- 
atively few gas particles. However, we feel that the 
level of convergence between the intermediate and 
high-resolution runs for objects with T > 1 keV is 
sufficient to give a good statistical representation 
of the net effect of feedback on the Lx-T relation. 

5. SUMMARY AND CONCLUSIONS 

In this paper we set out to extend the predic- 
tive power of current semi-analytic models (SAMs) 
of galaxy formation by investigating the effect of 
energy feedback from model galaxies on the prop- 
erties of intracluster gas. To achieve this objec- 
tive we have employed a novel hybrid technique in 
which a SAM is coupled to non-radiative hydro- 
dynamical simulations, thus guaranteeing that the 
source of feedback in our simulations is a realistic 
galaxy population. This is the first time such an 
approach has been adopted and is complementary 
to existing theoretical studies of galaxy groups and 
clusters based on self-consistent hydrodynamical 
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simulations. 

The main result to emerge from our work is that 
a large energy input from AGNs (on average, 35% 
of the available rest mass energy e r MBKC 2 ) is re- 
quired over the entire formation history of haloes 
in order to reproduce the observed Lx-T relation 
and hal o gas fractions. Th is supports the conclu- 
sion of iBower et al. ( 2008 ) derived using purely 
semi-analytic reasoning. 

We initially applied our method in a small 
L = 62.5k" 1 Mpc volume to explore how the bulk 
properties of groups and poor clusters are affected 
by different feedback components, concentrating 
on the Lx-T relation. The noteworthy results of 
this preliminary investigation are: 

1. The star formation and supernova feedback 
scheme employed in the Munich L-Galaxies 
SAM has a negligible effect on the entropy 
of the ICM, leading to an Lx-T relation 
that resembles the relation from a simula- 
tion with gravitational heating only. By 
contrast, the mean relation obtained from 
hydrodynamical simulations with radiative 
cooling and stellar feedback typically lies 
closer to the observed one, particularly on 
cluster scales. However, these simulations 
tend to over-produce stars, whereas the frac- 
tion of baryons locked-up in stars within the 
virial radius of our clusters is, on average, 
approximately 9%, in excellent agreement 
wth observational estimates. 

2. Incorporating a simple model for quasar 
mode AGN feedback in a simulation with 
L-Galaxies stellar feedback leads to X-ray lu- 
minosities considerably in excess of observed 
values, if the thermal coupling efficiency is 
ef = 0.05. The choice ef = 0.05 ensures 
the local black hole-bulge mass relation is 
recovered in hydrodynamical simulations in- 
cluding models for black hole growth and 



assoc iated quasar activity (|Di Matteo et al 



20051 ). However, our results suggest that 



simulations which e mploy the quasa r mode 



feedback scheme of iDi Matteo et al.l (|2008h 



would be unable to explain the observed 
scaling of X-ray luminosity with tempera- 
ture. 

We find that reproducing the desired steep- 
ening of the Lx-T relation on group scales 



actually requires a much larger coupling ef- 
ficiency: ef = 0.4. In this case the hot 
gas fractions of our simulated clusters also 
broadly agree with observations. This indi- 
cates that the balance between quasar and 
radio mode feedback needs to be adjusted 
in L-Galaxies if the model is to simultane- 
ously account for the observed properties of 
galaxies and the ICM. 



We ha ve implemented the recent lBower et al. 
(|2008h AGN feedback model in a hybrid 
simulation with stellar feedback from L- 
Galaxies. In this model, radio mode feed- 
back can eject X-ray emitting gas from cen- 
tral regions of a halo, reducing the gas den- 
sity and thus X-ray luminosity. This feed- 
back mechanism is more efficient in lower- 
mass systems, causing the Lx-T relation to 
become steeper on group scales as desired. 
Indeed, we find that the relation obtained 
from our simulation agrees well with obser- 
vational data on the mass scales probed. 
This is also true of the hot gas fractions of 
our simulated objects, demonstrating that 
significant AGN heating is a key ingredient 
in shaping groups and clusters. 



Once we had established a model capable of 
reproducing the observed Lx-T relation and gas 
fractions of groups and poor clusters, we investi- 
gated whether this best-fit model could also ex- 
plain the properties of richer clusters by perform- 
ing a hybrid simulation in a larger L = 125h~ 
Mpc volume. Our results can be summarised as 
follows: 

1. The observed Lx-T relation is successfully 
recovered on all mass scales, apart from the 
occasional object that has an anomalously 
high temperature due to a recent injection of 
energy from the central AGN. The variation 
in scatter along the relation also compares 
favourably with observations, although there 
appears to be a lack of objects at the upper- 
luminosity edge of the observed Lx-T rela- 
tion since we do not form CC systems in our 
non-radiative simulations. 

2. The gas fractions of our groups and clusters 
are in broad agreement with observational 
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data, displaying a similar degree of scat- 
ter. We find that AGN feedback significantly 
lowers the hot gas fraction in groups and 
poor clusters. This is because low-mass sys- 
tems have a shallow potential well and AGN 
heating can efficiently drive X-ray emitting 
gas from their central regions to their out- 
skirts. By contrast, massive clusters retain 
a greater fraction of their gas (up to 70% of 
the cosmic baryon fraction) since they have a 
much larger binding energy and AGN feed- 
back cannot effectively expel gas from the 
halo. 

The radial entropy profiles of our simulated 
clusters begin to flatten off at r ~ r 50 o , de- 
parting from the scaling S oc r observed in 
cluster outskirts. By contrast, observed en- 
tropy profiles typically flatten off at smaller 



radii r ~ 0-2r2nn (e.g iPonman et alJ l2003t 
Pratt et al ] 120061 ). Consequently, we tend 



to over-estimate the core entropy The core 
temperature is also higher than expected 
from observations of similar mass objects. 
However, profiles can only be reliably mea- 
sured for X-ray bright objects, which are 
probably a biased sample of the population 
of groups and poor clusters. Our model may 
provide a reasonable description of X-ray 
under-luminous systems for which is it not 
possible to extract profiles. Nevertheless, if 
we are to reproduce the observed profiles of 
the majority of X-ray bright groups and poor 
clusters, the core entropy and temperature 
must be reduced, potentially by including 
radiative cooling in our simulations. 

The work described in this paper represents 
merely the first stage in the development of our 
method. Even so, we have obtained several en- 
couraging results with this simple initial model. 
Following this success, we are currently resimulat- 
ing a sample of rich clusters (M500 ~ 10 15 /i _1 Mq) 
extracted from the Millennium volume with our 
best-fit feedback model. In such massive objects 
the gas cooling time is longer, so the lack of ra- 
diative cooling in our simulations is less of an is- 
sue. Indeed, preliminary results suggest that the 
density, temperature and entropy profiles we ob- 
tain are in better agreement with observed pro- 
files. These results will be presented elsewhere. 



Our primary goal in future work is to self- 
consistently incorporate radiative cooling into our 
hybrid approach, rather than relying on the simple 
cooling recipes employed in SAMs. These recipes 
usually assume that haloes have a spherically- 
symmetric isothermal gas distribution but, in gen- 
eral, neither of these assumptions will hold in 
hydrodynamical simulations. To circumvent this 
problem, we intend to fully couple SAMs to ra- 
diative simulations, so that the gas distribution in 
the simulation governs star formation, black hole 
growth and associated feedback in the SAM. This 
is a non-trivial task, requiring the simulation and 
SAM to be run simultaneously. Including gas cool- 
ing in our simulations should enable us to pro- 
duce more realistic ICM profiles in cluster cores, 
where the cooling time is short, since cooling acts 
to lower the gas entropy in such regions. 

Once our model is fully developed and tested, 
we shall conduct simulations in the full Millen- 
nium volume. The idea is to generate a large sam- 
ple of rich clusters that are consistent with the 
high-quality X-ray data available on these scales. 
An example of an important application of such a 
sample would be modellin g the selection fu nctions 
of X-ray surveys (e.g. ISahlen et al. 20091 ). This 
is essential to exploit the full power of clusters as 
cosmological probes of the expansion history of the 
universe. 

We thank V. Springel for supplying the merger 
tree software and G. De Lucia for providing the 
code for the L-Galaxies semi-analytic model. We 
are also grateful to A. Sanderson and A. Vikhlinin 
for making their observational data available to 
us. All simulations were performed using the 
Virgo Consortium Cosmology Machine at the In- 
stitute for Computational Cosmology, Durham. 
This work was supported by a Science and Tech- 
nology Facilities Council rolling grant. 
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Fig. 10. — Bolometric X-ray luminosity as a 
function of emission-weighted temperature for 
three hybrid simulations with our best-fit feedback 
model. To demonstrate the effect of resolution, the 
simulations are identical in all aspects except for 
the number of gas particles used, as detailed in the 
legend. 
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